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Abstract 

The  objective  of  the  present  study  is  to  investigate  the  application  of  the  baffle-blocked  flow  channels  for  enhancement  of  reactant  transport 
and  cell  performance  of  a  proton  exchange  membrane  fuel  cell  (PEMFC).  It  is  expected  that  due  to  the  blockage  effects  in  the  presence  of 
the  baffles,  more  fuel  gas  in  the  flow  channel  can  be  forced  into  the  gas  diffuser  layer  to  enhance  the  chemical  reactions  and  then  augment 
the  performance  of  the  PEMFC  systems.  The  effects  of  liquid  water  formation  on  the  reactant  gas  transport  are  taken  into  account  in  the 
present  modeling.  Predictions  show  that  the  local  transport  of  the  reactant  gas,  local  current  density  generation,  and  the  cell  performance  can 
be  enhanced  by  the  presence  of  the  baffles.  A  physical  interpretation  for  the  difference  in  the  baffle  effects  at  high  and  low  operating  voltages 
is  presented.  The  results  reveal  that,  at  low  voltage  conditions,  the  liquid  water  effect  is  especially  significant  and  should  be  considered  in  the 
modeling.  The  cell  performance  can  be  enhanced  at  a  higher  air  velocity  on  the  cathode  side,  by  which  the  cell  performance  can  be  enhanced 
and  occurrence  of  the  mass  transport  loss  can  be  delayed  with  the  limiting  current  density  raised  considerably. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  is  an  electro-chemical  system  that  can  realize 
the  direct  conversion  of  chemical  energy  of  reactants  (a  fuel 
with  an  oxidant)  to  electric  energy  with  high  efficiency  and 
high  environment  compatibility.  A  proton  exchange  mem¬ 
brane  fuel  cell  (PEMFC)  operates  at  significantly  low  tem¬ 
perature  than  other  types  of  fuel  cells  do.  It  is  a  potential 
candidate  of  the  power  source  in  the  near  future.  There  exist 
a  number  of  studies  on  the  transport  phenomena  and  the  per¬ 
formance  of  the  PEM  fuel  cells  in  the  open  literature.  Springer 
et  al.  [1]  proposed  an  isothermal,  one-dimensional  model  for 
the  proton  exchange  membrane  fuel  cells.  In  their  model, 
water  diffusion  coefficient,  electro- osmotic  drag  coefficient, 
water  sorption  isotherms,  and  the  membrane  conductivities 
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were  assumed  to  be  functions  of  membrane  water  content. 
Results  showed  that  the  membrane  resistance  can  be  signifi¬ 
cantly  raised  with  an  increase  in  current  density,  and  the  net 
water  flux  ratio  under  typical  condition  is  much  less  than  that 
under  the  fully  hydrated  membrane.  It  was  also  found  that  the 
resistance  then  reduces  for  a  membrane  of  thickness  small 
enough.  Kim  et  al.  [2]  presented  an  empirical  equation  to  fit 
the  experimental  data  of  cell  potential  and  current  density  for 
the  PEM  fuel  cells  at  various  operating  conditions  of  temper¬ 
atures,  pressures,  and  oxygen  compositions  in  the  cathode  gas 
mixture.  Singh  et  al.  [3]  employed  a  two-dimensional  model 
to  investigate  the  transport  phenomena  in  a  PEMFC.  In  their 
work,  the  attention  was  focused  on  the  transport  processes 
in  fuel  cells  in  order  to  improve  the  thermal  and  water  man¬ 
agement,  and  to  alleviate  the  mass  transport  limitations.  To 
examine  the  effects  of  temperature  and  gas  pressure  gradi¬ 
ents  on  fuel  cell  performance  and  water  management,  a  one¬ 
dimensional,  non-isothermal  model  was  proposed  by  Djilali 
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Nomenclature 

C  mass  fraction 

Cp  quadratic  drag  factor 

D  mass  diffusivity  (m2  s-1) 

F  Faraday  constant  (96,487  C  mol- 1 ) 

I  current  density  (A  m-2) 

J  transfer  current  density  (A  m-2) 

kp  permeability 

/c  gas  channel  length  (m) 

/g  gap  between  diffuser  layer  and  baffle  (m) 

N  number  of  baffles 

P  pressure  (Pa) 

R  universal  gas  constant  (8.314  J  (mol  K)-1) 

S  molar  production  rate  per  unit  area 

T  temperature  (K) 

u  velocity  in  the  v-direction  (ms-1) 

v  velocity  in  the  y-direction  (ms-1) 

V  operating  voltage  (V) 

w  oxygen  mass  fraction 

Zf  Spices  valence 

Greek  symbols 

a  transfer  coefficient  for  the  reaction 

s  porosity 

r]  surface  overpotential  (V) 

X  channel  width  fraction,  /g//c 

p  density  of  oxygen  (kg  m-3) 

o  ionic  conductivity  of  the  ionomer  (£2-1  m-1) 

r  tortuosity  of  diffusion  layer 

0  phase  potential 

v  viscosity  of  flow 

<fi  ionomer  phase  potential  (V) 

Subscripts 
a  anode 

c  cathode 

m  membrane 

d  diffuser  layer 

u  v-direct 

v  y-direct 

L  liquid  water 

eff  effective 


and  Lu  [4] .  It  was  found  that  the  temperature  distribution  in 
the  PEMFC  is  affected  by  the  water  phase  change  in  the  elec¬ 
trode.  The  peak  temperature  within  the  cell  occurs  at  lower 
operating  temperature  and  partially  humidified  reactants.  As 
a  result  of  enhanced  water  evaporation  with  thermal  gradients 
taken  into  account,  the  Darcy  isothermal  pressure  drop  can 
be  compensated  by  the  temperature  non-uniformity. 

In  the  study  of  the  PEM  fuel  cells,  the  thermal  and  water 
management  is  one  of  the  most  important  issues.  Nguyen  and 


White  [5]  proposed  a  model  of  the  water  and  heat  manage¬ 
ment  of  the  PEMFC  systems,  which  includes  the  effect  of 
electro-osmosis,  diffusion  of  water,  heat  transfer  from  solid 
phase  to  gas  phase,  and  latent  heat  as  water  evaporation  and 
condensation.  They  found  that  the  ohmic  loss  is  considerable 
at  high  current  density  and  the  voltage  loss  is  twice  amount 
of  that  of  the  cathode  electrode.  The  back-diffusion  of  the 
water  from  the  cathode  to  the  anode  is  not  enough  to  keep 
the  membrane  hydrated.  Therefore,  the  reactant  gas  at  the 
anode  needs  to  be  humidified.  In  their  work,  Voss  et  al.  [6] 
developed  a  relatively  novel  technique  of  water  management 
with  anode  water  removal  to  modify  the  water  concentra¬ 
tion  profile  or  gradient  of  the  proton  exchange  membrane 
to  augment  the  back-diffusion  rate  of  water  from  the  cath¬ 
ode  to  the  anode.  Therefore,  the  water  at  the  cathode  catalyst 
layer  diffuses  through  the  membrane  and  is  removed  via  the 
anode  reactant  gas  stream.  Yi  and  Nguyen  [7]  proposed  an 
along  the  channel  model  for  evaluating  the  effects  of  vari¬ 
ous  design  and  operating  parameters  on  the  performance  of 
a  PEMFC.  The  results  show  that  humidification  of  the  anode 
gas  is  required  to  enhance  the  conductivity  of  the  membrane, 
and  the  liquid  injection  and  higher  humidification  tempera¬ 
ture  can  improve  the  cell  performance  by  introducing  more 
water  into  the  anode.  Also  applying  higher  cathode  gas  pres¬ 
sure  helps  to  replenish  the  water  loss  by  electro-osmosis, 
thereby  making  the  membrane  more  conductive  and  thus  re¬ 
sulting  in  higher  cell  performance.  Baschuk  and  Li  [8]  devel¬ 
oped  a  mathematical  model  with  variable  degrees  of  water 
flooding  in  the  PEMFC.  Physical  and  electrochemical  pro¬ 
cesses  occurring  in  the  membrane  electrolyte,  the  cathode 
catalyst  layer,  the  electrode  backing  layer  and  the  flow  chan¬ 
nel  were  considered.  Compared  with  experimental  results, 
they  found  that  when  air  was  used  as  the  cathode  fuel,  the 
flooding  phenomena  are  similar  at  different  operating  condi¬ 
tions  of  the  pressures  and  temperatures.  When  cell  pressure  is 
increased  significantly,  the  water  flooding  in  the  electrode  be¬ 
comes  serious  and  leads  to  a  noticeable  reduction  in  the  power 
output. 

The  design  of  the  flow  field  in  bi-polar  plate  is  one  of  the 
crucial  factors  to  the  performance  of  a  PEMFC.  To  enhance 
the  liquid  water  transport  out  of  the  diffuser  layer,  a  new  flow 
channel  design  was  developed  by  Nguyen  [9].  This  design 
in  effect  has  converted  the  transport  of  reactant/product 
gases  to/from  the  catalyst  layers  from  diffusion  mechanism 
to  a  convection  mechanism.  Wood  et  al.  [10]  presented  the 
experimental  results  on  the  effectiveness  of  the  direct  liquid 
water  injection  scheme  and  interdigitated  flow  field  design 
for  proving  the  performance  of  PEMFC  s.  They  found  that 
the  interdigitated  flow  field  is  able  to  provide  higher  transport 
rates  of  reactant  and  products  to  and  from  the  inner  catalyst 
layers.  To  understand  the  effects  of  the  interdigitated  gas 
distributor  on  the  performance  of  these  fuel  cell  electrodes, 
Yi  and  Nguyen  [11]  proposed  a  two-dimensional  isothermal 
model  of  a  porous  electrode  to  simulate  the  hydrodynamics 
of  gas  flow  through  the  pore  volume  of  the  electrode  of  a 
PEMFC.  Based  on  the  predicted  results,  it  was  concluded 
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that  the  electrode  performance  is  improved  by  the  higher 
gas  flow  rate  through  the  electrode  improves.  Besides,  the 
average  current  density  decreases  with  an  increase  in  the 
electrode  thickness  or  shoulder  width  of  the  gas  distributor. 

He  et  al.  [12]  studied  the  effects  of  various  electrode  and 
flow  field  design  parameters  on  the  performance  of  the  cath¬ 
ode  of  a  PEMFC.  It  was  found  that  that  higher  differential 
pressure  between  inlet  and  outlet  channels  will  enhance  the 
electrode  performance.  Increasing  the  electrode  thickness  is 
equivalent  to  increasing  the  diameter  of  a  pipe  in  a  fluid  flow 
system  and  the  gas  bypassing  effect  becomes  more  signifi¬ 
cant.  Um  and  Wang  [13,14]  developed  a  multidimensional 
model  to  study  the  electrochemical  kinetics,  current  distri¬ 
bution,  fuel  and  oxidant  flow,  and  multicomponent  transport 
in  a  PEMFC  with  the  interdigitated  flow  field.  The  electro¬ 
chemical  behaviors  of  a  direct  methanol  fuel  cell  with  ser¬ 
pentine  flow  field  (SFF)  and  interdigitated  flow  field  (IFF) 
at  both  cathode  and  anode  were  investigated  by  Arico  et  al. 
[15].  The  experimental  results  showed  that  the  IFFs  signifi¬ 
cantly  enhance  the  mass  transport  inside  a  DMFC  and  allow 
achieving  the  higher  maximum  power  outputs  compared  to 
the  classical  serpentine  geometry.  The  larger  methanol  per¬ 
meation  through  the  electrolyte  determines  both  lower  volt¬ 
age  and  fuel  efficiencies  in  the  activation-controlled  region 
with  respect  to  the  SFF.  Recently,  Yan  et  al.  [16, 17]  presented 
the  studies  about  the  effects  of  flow  channel  designs  on  the 
steady  or  dynamic  gas  transport  and  cell  performance  of  PEM 
fuel  cells.  They  found  that  the  flow  distributor  geometry  has 
a  significant  influence  on  the  cell  performance.  In  addition, 
they  indicated  that  for  a  larger  channel  width  fraction,  a  faster 
dynamic  response  as  well  as  better  cell  performance  can  be 
attained. 

With  the  appropriate  design  of  the  flow  channel  of  the 
bipolar  plate,  the  efficiency  of  the  thermal  and  water  man¬ 
agement  in  PEMFC  systems  can  be  achieved.  Therefore  how 
to  design  an  optimal  flow  channel  in  the  bipolar  plate  is  one 
of  the  most  important  issues  in  the  study  of  the  PEMFC.  To 
install  a  tandem  array  of  baffles  in  fuel  flow  channels  is  a 
possible  strategy.  It  is  expected  that  due  to  the  blockage  ef¬ 
fects  in  the  presence  of  the  baffles,  more  fuel  gas  in  the  flow 
channel  can  be  forced  into  the  gas  diffuser  layer  to  enhance 
the  chemical  reactions  and,  in  turn,  augment  the  performance 
of  the  PEMFC  systems.  From  the  literatures  reviewed  above, 
it  is  obvious  that  the  performance  of  the  blocked  flow  chan¬ 
nels  with  baffles  installed  in  a  PEMFC  have  not  been  well 
examined  yet.  This  motivates  the  present  study.  The  objec¬ 
tive  of  the  present  study  is  to  investigate  the  performance  of 
the  blocked  flow  channel  with  baffles  in  a  proton  exchange 
membrane  fuel  cells  (PEMFC)  is  investigated  numerically. 
A  two-dimensional  numerical  model  is  developed  to  carry 
out  the  reactant  gas  transport  phenomena  and  cell  perfor¬ 
mance  with  baffle  effects  in  the  flow  channel  of  bipolar  plate. 
Effects  of  the  baffle  width  are  explored.  Additionally,  the  ef¬ 
fects  of  liquid  water  formation  on  the  reactant  gas  transport 
are  taken  into  account  in  the  modeling  and  examined  in  the 
analysis. 
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[T]  A  node  channel 
0  Anode  diffuser 
a  Anode  catalyst 
0  PEM 

[~5~|  Cathode  catalyst 
|~6~|  Cathode  diffuser 
[7]  Cathode  channel 
0  Bailie 


Fig.  1.  Two-dimensional  fuel  cell  model  of  a  PEMFC  system  with  baffle- 
blocked  flow  channel. 


2.  Analysis 

Fig.  1  is  a  schematic  diagram  of  the  two-dimensional  fuel 
cell  model  of  the  PEMFC  considered  in  the  present  work, 
which  consists  of  a  proton  exchange  membrane,  two  catalyst 
layers,  two  gas  diffuser  layers,  and  two  flow  channels.  To 
examine  the  effects  of  the  baffles  on  the  gas  transport  and  cell 
performance,  the  tandem  array  of  the  baffles  is  arranged  in  the 
cathode  flow  channel  in  order  to  force  the  fuel  gas  to  flow  into 
the  gas  diffuser  layer  and  catalyst  layer  at  the  cathode  side. 
For  simplification  of  the  analysis,  the  following  assumptions 
are  invoked:  (1)  the  gas  mixtures  are  considered  to  be  perfect 
gases;  (2)  the  gas  flow  is  steady  and  laminar  everywhere  in 
the  fuel  cell;  (3)  the  gas  diffuser  layers,  catalyst  layers,  and 
the  PEM  are  considered  to  be  isotropic  porous  media;  (4)  the 
heat  generated  by  the  chemical  reaction  is  neglected  in  this 
work. 

With  the  above  assumptions,  the  gas  transport  equations 
for  a  two-dimensional  PEMFC  system  can  be  depicted  as 
follows. 


Continuity  equation  : 


du  dv 
—  +  —  =0 
dx  d  y 


Momentum  equations: 


£eff 


du  du  \ 

u - 1-  v  — 

dx  dy  J 


£eff  aP 

p  dx 


+  F>£eff 
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Table  la 

Parameters  used  in  this  work 


Porosity  of  diffuser  layer 

0.4 

Porosity  of  catalyst  layer 

0.4 

Baffle  width  (m) 

0.000762 

Channel  length  (m) 

0.0762 

Channel  width  (m) 

0.000762 

Diffuser  layer  width  (m) 

0.000254 

Catalyst  layer  width  (m) 

0.0000287 

Membrane  width  (m) 

0.00023 

Operation  temperature  (K) 

353 

Operation  pressure  (atm) 

1 

Inlet  average  velocity,  vq  (ms-1) 

1.76 

(  dv  dv\  £eff  dP  ( d2v  d2v\ 

£eff  \uf. - H  V7T  I  =  —  b D£gff  (  + 

\  ox  dy  J  P  dy  \  dxz  dyz  J 

(3) 


£L  _  .  ^H20^c  8tffpRT°  (PU20  ~  Aat),  ^  Pu20  >  ^sat  ^ 

ke£effS(PSat  ~  Pu 2o)»  if  ^H20  <  ^sat 

where  M  is  the  molecular  weight  and  kc  and  kQ  are  the  con¬ 
densation  and  evaporation  rate  constants,  respectively.  The 
saturation  pressure  of  water  can  be  expressed  as 

psat  =  i0_2-1794+0-02953r_9-1837xl0_5:r2+L4454xl0_7  (6) 

In  addition,  the  saturated  rate,  s,  is  defined  as  the  ratio  of 
the  volume  of  pore  occupied  by  liquid  water  to  the  volume 
of  pore  in  the  porous  medium. 

For  this  reason,  the  effective  porosity  of  porous  media  is 
modified  to  account  the  liquid  water  effect: 

£eff  =  e(l  -  s)  (7) 


Species  concentration  equation  : 


£eff 


The  phase  potentials  in  the  catalyst  layers  satisfy  the  fol¬ 
lowing  equations: 


+ 


+  Sk  +  SL 


d 

dx 


at  anode 


In  the  momentum  equations,  Su  and  Sv  stand  for  the  source 
terms  based  on  the  Darcy’s  drag  forces  in  the  v  and y  directions 
imposed  by  the  pore  walls  on  the  fluid,  and  usually  cause  in  a 
significant  pressure  drop  across  the  porous  media.  The  details 
of  Su  and  Sv  for  different  layers  are  listed  in  Tables  la  and  lb. 

In  Eq.  (4),  the  effective  mass  diffusivity,  Deff,  can  be  writ¬ 
ten  as  Dk  Q ff  =  Dk£Ti  and  t*  is  the  tortuosity  of  the  pores  in 
porous  medium  and  Sk  is  the  production  rates  of  Mi  species 
in  gas  phase. 

To  comprehend  the  liquid  water  effect,  a  simplified  two- 
phase  model  is  used  to  describe  water  transport  in  the  PEM 
fuel  cell.  In  this  work,  we  account  for  the  liquid  water  effect  by 
modify  the  mass  diffusivity  due  to  the  liquid  water  filling  the 
pores  in  the  porous  media  and  the  liquid  water  generation  in 
the  species  concentration  equation.  When  the  partial  pressure 
of  water  vapor  is  greater  than  the  saturation  pressure  of  water 
vapor,  we  assume  that  water  vapor  may  condense  and  fill  the 
pore  in  the  porous  media.  The  source  term,  Sfi  in  the  species 
concentration  equation,  Eq.  (4),  representing  the  quality  of 
liquid  water  can  be  evaluated  [18]: 


d  /  d&\  d  (  d&\ 

—  <?m—  +  —  o-m—  =  —  7c  at  cathode 

dx  \  dx  J  dy  \  dy  J 


and  the  phase  current  density  satisfies  as  following  equations: 


3  ix  3 iy  .  , 

- - -  =  7a  at  anode 

A  1  A  J  Cl 

dx  dy 


(10) 


di 


X 


di 


H - -  =  jc  at  cathode 


dx  dy 


(ID 


The  governing  equation  of  the  proton  concentration  in  the 
membrane  can  be  written  as 


m 


BCk  ,  dCk 
u— - 1-  v 


dx 


dy 


-  D  I  ^  I  ^ 

ZF  (d2<P  d20 

+  ^0,effH+M/H+  +  Jy2 


(12) 


Table  lb 

Details  of  the  source  terms  in  the  governing  equations 


Gas  channel 
Gas  diffuser  layers 

Catalyst  layers 

Membrane 


J 


uz  +  vl  + 


k 

—  Z{Ch+F-  v<p- 

V 


du 

dx 


Sk 


Dk ,  eff 


Xi 


0 

7 

vs 

- v 


0 

0 


DkS 


n 


d,eff 


1.5 


vs 


£  „  nffCppv 


- v  — 

k 


P  v^p 


- — — V  Uz  +  vz 


vs 


slftCppv 


rs - V 

- V - _ — VM  +  Vz  + 

^p  \Ap 

kD  dv 

—  Zf CH+  F  ■  V<P  ■  — 
v  dx 


H2:  -  „  ,/a  02:  ——p—jc  H20: 


2  FC* 


AFCr 


2  FCc 


jc  Dks 


ZF  ^  ^  (d2<P  t  d20 

+  CH+  +  ~kyl 


r  i 


c,eff 


DkS 


X; 


m,eff 
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and  the  membrane  phase  potential  is  governed  by  the  follow¬ 
ing  equation: 


d 

dx 


(13) 


In  the  above  formulations,  Eqs.  (1)— (4)  and  (8)— (1 3)  form 
a  complete  set  of  governing  equations.  The  detailed  boundary 
conditions  are  similar  to  those  of  Soong  et  al.  [19].  Due  to 
the  space  limit,  they  are  not  presented  here. 


3.  Numerical  method 

Since  it  is  impossible  to  obtain  an  analytic  solution  of  the 
complex  convection-diffusion  problem  like  this  one,  it  will 
be  solved  by  a  finite  volume  method  using  a  collocated  cell- 
centered  variable  arrangement.  The  governing  equations  can 
be  expressed  in  the  form  of  a  generalized  transport  equation 

v  •  ( Pu(f>  -  r0v0)  =  (14) 

where  0  is  the  general  dependent  variable,  r $  the  exchange 
coefficient,  the  source  term,  u  velocity  vector,  and  p  the 
density.  With  the  discretization  of  the  governing  equations, 
the  coupled  finite-difference  equations  become 

<3P0P  =  «E0E  +  «W0W  +  ^N0N  +  <2S0S  +  S 0  (15) 

where  0p  is  the  value  of  0  at  the  current  point  P,  0e.  .  .0s 
stand  for  the  values  of  the  grid  points  adjacent  to  the  point  P, 
and  ap. .  .as  are  known  as  the  link  coefficients. 

In  this  work,  the  non-uniformly  distributed  grid  systems 
of  101  and  183  points  are  employed  in  the  x  and  y  directions 
(or  expressed  as  101  x  1 83),  respectively.  To  obtain  better  ac¬ 
curacy  in  the  numerical  computations,  grid  independence  is 
examined  in  the  preliminary  test  runs.  It  is  found  that  the  de¬ 
viations  of  the  local  current  density  between  the  predictions 
on  the  grids  of  101  x  183  and  201  x  285  are  less  than  2%.  Ac¬ 
cording,  the  computations  on  the  grid  system  of  101  x  183 
points  seems  sufficient  to  understand  the  behaviors  of  the 
reactant  gas  transport  in  a  PEMFC.  To  further  check  the  ad¬ 
equacy  of  the  numerical  scheme,  at  first,  the  results  of  the 
two-dimensional  results  without  the  baffle  effect  are  evalu¬ 
ated,  as  shown  in  Fig.  2.  It  is  clearly  seen  from  Fig.  2  that 
the  present  predictions  agree  reasonably  with  the  experimen¬ 
tal  data  of  Chen  [20].  The  above  preliminary  runs  confirm 
that  the  present  model  and  the  numerical  method  used  are 
generally  appropriate  in  analysis  of  the  present  problem. 


4.  Results  and  discussion 

With  the  tandem  arrangement  of  the  baffles  in  the  chan¬ 
nel  of  the  cathode  side,  the  fuel  gas  will  be  forced  to  flow 
into  the  gas  diffuser  layer  and  catalyst  layer  to  enhance  the 
chance  of  chemical  reaction  at  the  catalyst  layer  as  well  as 
the  performance  of  the  fuel  cell  system.  To  examine  the  baffle 
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Fig.  2.  Comparison  of  the  predicted  I-V  curve  without  baffle  effect  and  the 
experimental  data  of  Chen  [20] . 

effects  on  the  cell  performance,  the  polarization  (I-V)  curves 
with  various  baffle  width  and  baffle  numbers  are  presented 
in  Fig.  3(a)  and  (b),  respectively.  For  comparison,  the  results 
without  liquid  water  effects  in  the  model  are  also  presented. 
An  overall  inspection  of  Fig.  3  indicates  that  at  the  conditions 
of  the  high  operating  voltages,  the  baffle  effects  on  the  overall 
cell  performance  are  negligibly  small  as  compared  with  those 
without  baffle  effect.  At  low  voltage  conditions,  however,  the 
baffle  effects  on  the  I-V  curves  become  important.  In  addi¬ 
tion,  the  cell  performance  is  enhanced  with  increase  in  the 


Fig.  3.  Effects  of  baffle  width  and  baffle  numbers  on  the  cell  performance: 
(a)  baffle  width  effects  and  (b)  baffle  number  effects. 
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Fig.  4.  Effects  of  baffle  width  on  the  local  current  density  distributions:  (a) 
V=  0.7  V  and  (b)  V=0.2V. 


Fig.  5.  Effects  of  baffle  number  on  the  local  current  density  distributions: 
(a)  V =0.7  V  and  (b)  V  =  0.2  V. 


baffle  width  and/or  baffle  numbers.  It  is  clearly  seen  that  the 
higher  limiting  current  density  is  noted  for  a  PEMFC  with  a 
greater  baffle  width  and/or  baffle  numbers.  It  is  also  disclosed 
that,  at  high  operating  voltage  conditions,  the  deviations  in 
the  cell  performance  between  the  results  with  and  without 
consideration  of  liquid  water  effects  are  small.  This  means 
that  the  fuel  transport  in  the  PEMFC  can  be  safely  treated 
as  gas  phase  only  at  a  high  voltage  condition.  Whereas,  at 
low  voltage  conditions,  the  liquid  water  effects  on  the  cell 
performance  are  remarkable  and  cannot  be  neglected  in  the 
modeling.  This  confirms  the  fact  that  the  mass  transports  at 
low  V  are  significant  and,  in  turn,  more  water  is  generated  in 
the  catalyst  layer  of  the  cathode  side.  Therefore,  two-phase 
flow  effects  should  be  considered  in  the  low  operating  voltage 
conditions. 

The  above  evidences  demonstrated  that  the  appearance 
of  the  baffles  in  the  flow  channel  has  a  considerable  impact 
on  the  cell  performance.  It  is  very  significant  to  examine 
the  baffle  effects  on  the  local  transport  characteristics  of  the 
PEMFC,  especially  at  the  operating  conditions  of  low  volt¬ 
age.  To  reach  this  end,  Figs.  4  and  5  show  the  effects  of  baffle 
width  and  baffle  numbers  on  the  local  current  density  distri¬ 
butions,  respectively.  In  these  plots,  the  baffles  are  located  at 
the  positions  of  i/(N+  l)th  length  from  the  cathode  inlet.  For 
example,  for  case  of  one  baffle  (i.e.,  N=  1),  the  baffle  location 


is  on  the  location  of  the  l/2th  length  from  the  inlet,  that  is, 
at  the  central  position  of  the  channel  length.  For  N  =  3,  the 
three  baffles  are  located  at  the  positions  of  1/4,  1/2  and  3/4th 
length  from  the  inlet,  respectively.  It  is  clearly  observed  from 
Fig.  4  that  except  the  peak  regions,  the  local  current  density 
decreases  with  the  axial  location  Y.  However,  around  the  baf¬ 
fle  locations,  there  exist  peaks  in  the  distributions  of  the  local 
current  density.  It  means  that  better  cell  performance  is  gen¬ 
erated  around  the  baffle  locations.  This  can  be  made  plausible 
by  noting  the  fact  that  as  the  fuel  gas  is  blocked  by  the  baffle 
in  the  channel,  more  fuel  gas  will  be  forced  to  turn  into  and 
penetrate  through  the  gas  diffuser  layer  and  catalyst  layer, 
which  in  turn,  causes  a  significant  enhancement  of  chemical 
reaction  at  the  catalyst  interface.  The  plots  show  local  current 
density  distributions  of  a  periodic  pattern  with  peak  value  ap¬ 
pearing  in  the  regions  around  the  baffles  and  reducing  along 
the  channel.  The  best  enhancement  does  always  occur  at  the 
region  around  the  first  baffle.  A  careful  examination  of  the 
cases  with  V  =  0.7  V  in  Fig.  4(a)  discloses  that,  the  local  cur¬ 
rent  density  monotonically  increases  with  an  increase  in  the 
baffle  width  in  the  upstream  region  quite  near  the  entrance; 
while,  at  the  downstream  region  (Y >  0.065  m)  of  the  blocked 
channel  with  wider  baffles,  remarkable  reductions  in  current 
density  are  found  and  the  reduction  becomes  more  serious  as 
the  baffle  number  N  increases.  The  dramatic  degradation  in 
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local  current  density  in  the  downstream  region  of  a  blocked 
channel  is  a  reflection  of  the  more  efficient  fuel  transport  and 
the  chemical  reaction  in  the  upstream.  Consequently,  these 
two  counter-effects  render  the  resultant  current  density  or 
performance  of  the  fuel  cell  insensitive  to  the  application  of 
the  baffles.  For  the  cases  with  V=0.2  V  in  Fig.  4(b),  abrupt 
increases  in  current  density  appear  due  to  the  presence  of  the 
baffles,  and  the  high-/ regions  beneath  the  baffles  expend  with 
an  increase  in  the  baffle  width.  Although  the  values  of  the  cur¬ 
rent  density  in  the  regions  between  the  neighboring  baffles 
reduce  with  the  increasing  baffle  width,  in  sum,  noticeable 
enhancement  in  the  cell  performance  can  be  still  retained.  The 
global  effects  addressed  above  for  the  cases  with  the  baffles 
of  various  width  at  the  operating  conditions  of  0.7  and  0.2  V 
can  be  validated  in  the  I-V  curves  in  Fig.  3(a).  Besides,  the 
predictions  reveal  that  consideration  of  the  liquid  water  in  the 
simulation  leads  to  the  relatively  lower  current  density  than 
those  in  the  cases  without  liquid  water  effect. 

Fig.  5  presents  the  effects  of  baffle  numbers  on  the  local 
current  density  distributions  for  the  baffle  width  W=  4  and  at 
the  operating  voltages  of  0.7  and  0.2  V.  Generally,  the  quali¬ 
tative  trend  of  the  data  in  Fig.  5  can  be  interpreted  in  a  similar 
way  as  that  we  did  for  the  cases  of  various  baffle  width  ef¬ 
fects  in  Fig.  4.  The  major  difference  between  two  is  that  the 
enhancement  in  local  current  density  in  Fig.  4  is  achieved  by 
enlarging  the  extend  of  the  high-/  regions  via  an  increase  in 
the  baffle  width,  whereas  in  the  cases  shown  in  Fig.  5,  the  local 
current  density  or  performance  is  enhanced  by  increasing  the 
number  of  the  high-/ regions.  The  effects  of  the  baffle  number 
effects  at  the  operating  condition  of  low  voltage  (V=0.2  V) 
are  shown  in  Fig.  5(a).  Near  the  entrance,  the  effects  of  baffle 
on  the  local  current  density  are  similar  to  those  at  the  high 
operating  voltage.  For  the  downstream  region,  however,  the 
peak- value  regions  in  the  local  current  density  plots  become 
narrow  and  shortened.  This  can  be  attributed  to  the  more  pro¬ 
nounced  electro-chemical  reaction  at  the  operating  condition 
of  low  voltage.  Therefore,  oxygen  concentration  at  the  cath¬ 
ode  side  would  be  reduced  considerably,  especially  for  the 
cases  with  a  greater  number  of  baffles.  Accordingly,  in  the 
downstream  region  of  the  blocked  channels,  the  high-/ region 
around  the  baffle  becomes  narrow. 

The  distributions  of  oxygen  mass  flow  rate  along  the  in¬ 
terface  between  the  gas  diffuser  layer  (GDL)  and  catalyst 
layer  (CL)  with  various  arrangements  of  the  baffles  are  pre¬ 
sented  in  Fig.  6.  Since  the  power  output  is  the  consequence  of 
the  electro-chemical  reaction,  the  consumption  of  the  oxygen 
along  the  GDL-CL  interface  can  be  considered  as  an  index 
of  the  cell  performance.  A  higher  mass  flow  rate  indicates 
a  higher  current  density,  i.e.,  a  better  cell  performance.  It  is 
apparent  that  the  distributions  of  the  oxygen  mass  flow  rate 
are  similar  to  those  of  the  local  current  density.  This  confirms 
the  adequacy  of  the  present  study. 

In  the  study  of  reactant  gas  transport  in  the  PEMFC,  un¬ 
derstanding  of  the  detailed  local  distributions  of  the  oxygen 
at  the  cathode  side  along  the  GDL-CL  interface  is  important 
to  the  design  of  a  PEMFC.  To  this  end,  the  effects  of  baffle 


Fig.  6.  Effects  of  baffle  width  and  baffle  number  on  the  local  distributions 
of  oxygen  mass  flow  rate. 

width  and  baffle  numbers  on  the  local  distributions  of  oxygen 
mass  fraction  along  the  GDL-CL  interface  for  V  =  0.2  V  are 
presented  in  Fig.  7,  in  which  it  is  observed  that  the  presence 
of  the  baffles  in  the  flow  channel  has  a  considerable  impact  on 
the  local  distributions  of  oxygen.  Around  the  baffle  regions, 
the  oxygen  distributions  show  a  peaked  increase,  which  is 
a  consequence  of  the  channel  blockage  and  its  effect  on  the 
enhancement  of  fuel  gas  transport  we  have  addressed  previ¬ 
ously.  Therefore,  the  oxygen  distributions  show  peaked  vari¬ 
ations  around  the  baffles.  In  Fig.  7(a),  for  the  peaked  variation 
around  the  first  baffle,  the  influenced  regions  due  to  the  pres¬ 
ence  of  the  baffle  increases  with  the  increased  baffle  width. 
But  for  the  third  peaked  variation,  the  smallest  influenced  re¬ 
gion  and  peak  value  in  the  oxygen  concentration  are  noted 
in  the  case  with  a  widest  baffle,  in  which  the  chemical  reac¬ 
tion  along  the  catalyst  in  stronger  near  the  entrance.  But  as 
the  fuel  gas  goes  downstream,  the  oxygen  concentration  be¬ 
comes  lower  due  to  the  faster  consumption.  Therefore,  in  the 
downstream  region,  the  oxygen  concentration  may  become 
lower  around  the  baffle  for  the  case  with  a  wider  baffle.  As  to 
the  baffle  numbers  effects  on  the  local  oxygen  distributions, 
Fig.  7(b)  presents  three  typical  predictions  of  N=  1,  3,  and  5. 
Generally  speaking,  as  those  presented  in  Fig.  7(a),  the  pres¬ 
ence  of  the  baffles  leads  to  abrupt  and  high-peak  variations  of 
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Fig.  7.  Effects  of  baffle  width  and  baffle  numbers  on  the  local  distributions 
of  oxygen  mass  fraction  along  the  interface  between  the  GDL  and  catalyst 
layer  for  V  =  0.2  V. 

the  local  oxygen  mass  fraction  in  the  regions  corresponding 
to  the  locations  of  the  baffles.  Comparing  the  local  oxygen 
mass  fraction  in  the  region  beneath  the  baffle,  i.e.,  around 
the  1/2  channel  length,  in  the  three  cases,  it  is  shown  that 
the  local  oxygen  mass  fraction  reduces  with  increasing  baffle 
number. 

In  analysis  of  the  PEMFC,  the  airflow  rate  on  the  cathode 
side  is  one  of  the  key  operating  parameters  to  the  cell  per¬ 
formance.  To  examine  this  effect,  Fig.  8  shows  the  effects  of 


Fig.  8.  Effects  of  air  velocity  on  the  cathode  side  on  the  cell  performance. 


air  velocity  on  the  cathode  side  on  the  performance  curves 
of  the  PEMFC.  In  this  work,  the  reactant  inlet  velocity  of 
the  values  0.44,  0.88  and  1.76  m  s-1  are  considered.  The  pre¬ 
dictions  in  Fig.  8  disclose  that  a  better  cell  performance  is 
found  for  a  case  with  a  higher  air  velocity.  The  high  velocity 
implies  sufficient  supply  of  the  reactant  gas,  which  reduces 
the  mass  transfer  overpotential,  delay  the  occurrence  of  the 
mass  transport  loss  and  then  increases  the  limiting  current 
density  considerably.  Besides,  the  model  without  considera¬ 
tion  of  liquid  water  effect  over-predicts  the  simulated  results 
of  the  I-V  curves. 


5.  Concluding  remarks 

With  an  appropriate  design  of  the  flow  channel  of  the  bipo¬ 
lar  plate,  the  thermal  and  water  management  of  the  PEMFC 
can  be  achieved  efficiently.  To  reach  this  end,  a  new  con¬ 
cept  of  the  flow  channel  design  with  baffles  is  proposed  and 
analyzed  in  this  work.  To  validate  the  concept,  the  detailed 
gas  transport  phenomena  and  cell  performance  with  baffle  ef¬ 
fects  in  the  flow  channel  are  numerically  investigated.  What 
follows  are  the  major  findings. 

1 .  B  oth  the  reactant  transport  and  cell  performance  can  be  en¬ 
hanced  by  the  presence  of  the  baffles  in  the  flow  channel  of 
the  bipolar  plate,  especially  at  the  operating  conditions  of 
low  voltage.  The  beneficial  baffle  effects  become  increas¬ 
ingly  remarkable  with  increasing  width  and/or  number  of 
baffles  in  the  tandem  array. 

2.  Predictions  demonstrate  the  baffle  effects  on  the  cell 
performance  are  dependent  of  operating  conditions.  At 
high  operating  voltages,  the  baffle  effects  enhance  the 
gas  fuel  transport  and  raise  the  local  current  density 
in  the  upstream  region  but  reduce  them  in  the  down¬ 
stream  part  of  the  flow  channel.  The  resultant  influ¬ 
ence  of  the  presence  of  the  baffles  becomes  insignif¬ 
icant  due  to  offset  of  the  two  counter-effects.  At  low 
operating  voltages,  however,  abrupt  increases  in  local 
current  density  appear  at  the  locations  beneath  the  baf¬ 
fles.  Although  the  values  of  the  current  density  in  the 
regions  between  the  neighboring  baffles  reduce,  notice¬ 
able  enhancement  in  the  resultant  cell  performance  is  still 
retained. 

3.  The  predictions  obtained  by  using  the  models  with  and 
without  liquid  water  effects  are  of  deviations,  especially 
at  the  operating  conditions  of  low  voltage.  Under  this  sit¬ 
uation,  the  liquid  water  effect  is  significant  and  should  be 
taken  into  account  in  the  modeling. 

4.  With  a  higher  air  velocity  on  the  cathode  side,  the  lo¬ 
cal  current  density  is  increased  and  the  cell  perfor¬ 
mance  can  thus  be  enhanced.  The  high  velocity  of  the 
reactant  gas  can  reduce  the  mass  transfer  overpoten¬ 
tial,  delay  the  occurrence  of  the  mass  transport  loss, 
and  then  increase  the  limiting  current  density  conside¬ 
rably. 
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